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A detailed theoretical interpretation of the Josephson interference experiment between Sr2Ru04 
and Pb reported by Kidwingira et al [J is given. Assuming chiral p-wave pairing symmetry a 
Ginzburg-Landau theory is derived in order to investigate the structure of domain walls between 
chiral domains. It turns out that anisotropy effects of the Fermi surface and the orientation of the 
domain walls are essential for their internal structure. Introducing a simple model for a Josephson 
junction the effect of domain walls intersecting the interface between Sr2Ru04 and Pb is discussed. 
It is shown that characteristic deviations of the Fraunhofer interference pattern for the critical 
Josephson current as a function of the magnetic field occurs in qualitative agreement with the 
experimental finding. Moreover the model is able also to account for peculiar hysteresis effects 
observed in the experiment. 

PACS numbers: 74.20.-z 71.18.+y 



I. INTRODUCTION 

The aim to identify the symmetry of the supercon- 
ducting phase in Sr2Ru04 has stimulated numerous ex- 
periments since more than a decade 0, Q . Most of them 
provide strong evidence that the pairing state has the so- 
caUed chiral p-wave symmetry. This state is the analog 
of the A-phase of superfluid ^He, which is a spin triplet 
state breaking time reversal symmetry with an angular 
moment along the fourfold c-axis of the tetragonal crystal 
lattice of Sr2Ru04. The full gap function is a 2 x 2-matrix 
in spin space, Aj, and can be represented by the d-vector 
in the case of spin-triplet pairing: 

d(k) = AQz{kx i iky) — ~tr ^Aj:iay&^ , (1) 

with Aq as the gap magnitude. SQUID interferome- 
ter experiments probing the internal phase structure of 
the Cooper pairs are consistent with odd-parity (p-wave) 
pairing [J]. NMR Knight shift measurements are com- 
patible with equal-spin spin triplet pairing with the spin 
axis lying in the basal plane of the tetragonal crystal 
lattice, i.e. d || z Q. Muon spin relaxation studies 
show enhanced internal magnetism in the superconduct- 
ing phase suggesting a state with broken time reversal 
symmetry [6l|. Similarly, recent Kerr effect studies im- 
ply the presence of an orbital magnetic moment pointing 
along the c-axis Q- On the other hand, the search for 
the spontaneous magnetization at the surface of samples, 
as expected for a chiral p-wave state, by highly sensitive 
scanning probes has only given negative results so far 

A further extraordinary property of the chiral p-wave 
state is its two-fold degeneracy which can lead to do- 
mains, which are distinguished by the orientation of the 
orbital angular momentum of the Cooper pairs, parallel 
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or antiparallel to the c-axis. Early on the question of the 
formation of domains of the two states (angular momen- 
tum up and down) has been discussed, but no experi- 
mental indications of domains had been reported, until 
recently. Kidwingira et al investigated carefully the inter- 
ference pattern of a Josephson junction between Sr2Ru04 
and the conventional superconductor Pb in a magnetic 
field T|. They interpreted these patterns as the result of 
domain walls intersecting the extended Josephson junc- 
tion and so giving rise to a spatial variation of the Joseph- 
son phase along the junction. It has been speculated 
in the past that domain walls intersecting the interface 
between a chiral p-wave and conventional superconduc- 
tor could influence the Josephson effect, as the Joseph- 
son current-phase relation could be different for the two 
types of domains ■ Such a property could lead to an 
intrinsically inhomogeneous junction and would alter the 
interference pattern. This situation bears some similar- 
ity with the 45° asymmetric interface in oriented films of 
high-Tc superconductors where faceting of the boundary 
introduces a random switch of Josephson phase by ±7r 
along the interface [llj . This type of boundaries display 
unusual interference patterns of the critical current in a 
magnetic field and generate spontaneous flux pattern on 
the interface. Unlike in the standard case the maximum 
of the critical current is here usually not located at zero 
field. 

This finding could eventually provide indirect evidence 
for the presence of domain walls in the p-wave supercon- 
ductor. The discussion of this problem requires a detailed 
analysis of the domain wall structure. For interpretation 
of their results Kidwingira et al introduced several types 
of domain walls which could give rise to strong variations 
of the Josephson phase fl*!. While they are possible on a 
topological level, not all the proposed domain walls are 
energetically stable. Indeed in the most simple approach 
assuming full rotation symmetry for the Fermi surface 
around the z-axis of the p-wave superconductor none 
of the stable domain walls give rise to any shift of the 
Josephson phase and the presence of domain walls could 
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go basically unnoticed in the interference experiment . 
However, as we will show below a more careful analysis 
taking more general conditions of the electronic spectrum 
into account leads to stable and metastable domain walls 
which can give rise to non-trivial intrinsic phase patterns 
in a Josephson junction and reproduce the experimen- 
tally observed anomalies in the interference pattern. 



II. 



GINZBURG-LANDAU FREE ENERGY 



Our study on the Josephson effect involves both a con- 
ventional s-wave and a chiral p-wave superconductor, 
representing the experimental arrangement of Pb cou- 
pled to Sr2Ru04 T]. We first introduce here the ba- 
sic order parameters and their corresponding Ginzburg- 
Landau theories which will be used later to discuss the 
structure of domain walls in the chiral p-wave state and 
their influence on the Josephson effect in configurations, 
as shown in Fig. [1] 

The conventional superconductor is described by a 
scalar order parameter 'il'{r) for the spin-singlet pair- 
ing state of highest possible symmetry (" s-wave" pairing 
state). The Ginzburg-Landau free energy functional has 
the standard form. 



(2) 

with = V — i^A as the gauge-invariant derivative 
and 7 = 2e/fi,c = 27r/$o is the vector potential with 
the magnetic field S = V x A), $0 is the magnetic flux 
quantum, as{T) — a'^{T — Tcs), bg and Kg are parameters 



qua i 

m 



The chiral p-wave phase requires a two-component or- 
der parameter t/ = {r]x,r]y) with 



d{k) = z{r]xkx + rjyky) 



(3) 



which belongs to the irreducible representation Eu of the 
tetragonal point group of Sr2Ru04. The free energy 
functional has then the following general form 



Tp[r]^,r]y,A] = / 
Jv 



ap{T)\rj\^ + bi\v\ 



satisfy the relations: 62 > 0, 62 > &3 and 46i > 62 — 63- 
Minimizing the free energy functional with respect to I770I 
in the homogeneous case, we obtain the uniform phase^ 



T7± = ?7o(l,±«) , \vo\" 



4bi -62 + 63 
such that the corresponding gap function is 

d±{k) = r]o{T)z{kx ± iky) . 



(4) 



(5) 



The two states with opposite relative sign are degener- 
ate and violate time reversal symmetry, as the operation 
of time reversal K yields Kd± = d'^. = dip. These two 
states can form domains. In the following we will first an- 
alyze the structure of domain walls between such domains 
of opposite chirality and their influence on the Josephson 
effect in geometries as shown in Fig IT] 
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FIG. 1: Josephson Junction between the s-wave state and the 
chiral p-wave state with a domain wall parallel to the crystal 
a;-axis between the two degenerate states: t/^ oc (±l,i). The 
dimensions of the junction are Aa; x Ay x Ajz = 2Lx x Lx L^. 



III. DOMAIN WALL STRUCTURE 



+Xi(|i?,?7,|2 + \Dy7^y\'')+K2{\D,riy\'' + \DyT^^\') 

+ {Kz{D,ri,)*{DyT^y) + K4{D,f]y)*{DyT],) + c.c.} 



+K5{\D,r,,\^ + \D,r]y\ 



?1 

Stt 



The coefhcients bi and Ki are material-dependent param- 
eters and ttp = a'p{T — T^p) [13|]. In order to stabilize the 
chiral p-wave state in the bulk the coefficients have to 



Before addressing the influence of the domain wall on 
the Josephson effect it is necessary to analyze here the 
structure of a domain wall in detail. Our aim is to show 
which types of domain walls are energetically favorable. 
It will turn out that the anisotropy of the Fermi surface 
plays an important role and with this also the orientation 
of the domain wall. For a convenient discussion it will be 



^ We note that the corresponding bulk energy density is given by 
fhomoglvo] = (-4fei + b2 - b-i)rj^ = ap-ql = -HI/Stt, which 
defines the critical magnetic field He- 
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advantageous to change the order parameter representa- 
tion of the p-wave superconductor and to formulate the 
Ginzburg-Landau functional also in rotated coordinate 
frames. 



In this way we can use the angle to define the coordi- 
nate frame. In the following we will omit the primes and 
always assume that we describe the domain wall in the 
corresponding frame. 



A. Free energy functional formulation 

In a first state we express now the free energy func- 
tional in a new form which simplifies the rotation of the 
coordinate frame. We then always define the domain 
wall as the ^-x-plane and the spatial variation of the 
order parameter occurs only along the corresponding 
y-axis perpendicular to this plane. For simplicity we 
keep the z-axis as fixed along the crystaline z-axis. 



We introduce the parametrization |14l |. 

V± = ^(=^^2; - illy) , D^± iDy , 

which inserted into the free energy leads to 



(6) 



2ap(|r;+p + |r;_n 



In a weak-coupling approach the parameter v = {{v*) — 
+ i^x'^'y)) is a measure for the anisotropy 
of the Fermi surface (i^ = for a cylindrical Fermi sur- 
face and |:^| = 1 for a square-shaped Fermi surface) where 
Vx,y,z are the components of the Fermi velocity and (•) 
defines the average over the Fermi surface. The coeffi- 
cients in Eq.Q satisfy the following relations 



bi = , 
6, = -felM^ 



K3 



. (7) 



where b and K are material dependent parameters. 

This form of the free energy functional allows us now 
to deal easily with the rotation of the reference frame 
around the z-axis: = (cccos^ -I- t/sin0,2/cos0 — 

xsmO) with the angle 9 relative to the original x-axis 
in the tetragonal crystal. For the new coordinate frame, 
the order parameter and the gradients are transformed 
as ri± = e^*^77j_ and D± = e'^^^D'^. When we express 
([7]) in the new coordinates, we have only to modify the 
two terms containing the parameter v by phase factors, 

bi^ e T]\_ ri_ + c.c. , 



K 



B. Variational ansatz for the domain wall structure 

We now turn to the structure of the domain wall which 
will depend qualitatively on the choice of parameters in 
the free energy. We choose here a variational approach to 
discuss behavior of the order parameter around the do- 
main wall with the following ansatz which is most useful 
to eventually obtain the key information relevant for the 
Josephson effect. 



= Va^^^* cosx and 77- = 7706*'^ 



smx, 



(8) 



with the boundary conditions within a given reference 
frame, 







-00 
-00 



(9) 



where we restrict the spatial dependence of the order pa- 
rameter to the function xiu) and use the relative phase 
a — — (/)- as a further (variational) parameter. In this 
way the domain wall appears as an interface between the 
two superconducting domains which is additionally char- 
acterized by a phase difference a, similar to a Josephson 
junction or weak link [T^ . 

For the later discussion of the Josephson effect with an 
s-wave superconductor, it will be useful to return to the 
order parameter components (rj^, rjy) which are expressed 
as 



e'"^+ cos X 



cosx 



sinx) , 
- sinx) • 



(10) 



The chosen boundary conditions in © correspond to the 
situation, 



v{y = -00) = ?7o(-l,*)c**" = V 
q{y = +00) = 77o(+l,^)e'^+ = 77+ 



(11) 



for which, in the parametrization {VxtVv) — 
(|?7£c|e*'^^, Ijyyle"^"), the phase shifts of the order pa- 
rameters are given by 



^4>x = (+00) - 01 (-00) = a - TT 

/S.<f)y = (j)y{ + 00) — (j)y{ — = «. 



(12) 



Here a plays the role of the total phase difference of 
the order parameter and also determines the current flow 
through the domain wall analogous to a Josephson junc- 
tion. 
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With this variational ansatz we rewrite the free energy, spending three equations: 



+ ^(1 + i^cos(2a - 46'))sin2(2x) 



(V X Af 



i4r77o||-DcosxP + l-Dsinxl 



i(e-'"p+ cosx)*P- sinx) 



+i/e*("-4e)(£)^ sinx)*(I?- cosx) + c.c.)| 



where we substitute Op = — using (jH) and ([7]). As- 
suming homogeneity along the z- and x-axis, we take 
X = x(y) a^nd A = ^(y) with — 0, the free energy 
can be written as 



- bvl + Kr^l\^^ sin' 2x + {dyx? 



WiAl + AD + C+ sin2x (7^(4 - - ^^vX) 
+2S--i'A^Ay sin2x + 279j,x {S+Ay - C_yl^) 
, {dyA^f 



Stt 



where — K/2briQ defines the coherence length. Addi- 
tionally we introduced 



Q = 1 + i^cos(2a- 40), 

C± = (cos a ± I' cos(q - 46'))/2, 

S± = (sin a±i' sm{a - 4:9))/ 2 . 



(13) 



which are the only coefficients depending on the 
anisotropy parameter v, the phase a and the angle 9 
of the domain wall orientation relative to the crystalline 
main axis (crystal a;-axis). 

Certain symmetries become immediately obvious here. 
The free energy is invariant under the rotation 9 — > 
9zLtt/2 and the operation {9, v) — >■ (6'±7r/4, —v) leaves the 
coefficients unchanged. It is also clear that the phase a 
is closely linked to the orientation of the domain wall. A 
further aspect of symmetry is connected with the opera- 
tion (0, a) {~9, ~a), leading to Q +Q , C± +C± 
and S± — >■ ^S±. 

The variational minimization of the free energy func- 
tional with respect to Xj A^ and Ay leads to the corre- 



- cos 2x 



a 



C-dyA,) + C+sm2xd'yX 



-2S-j'A^Ay 



i^HAl-Al)-{dyx)')+2S-j 



'«^Co7 - 7(1 - C+ sin2x)A^ = S^-/Ay sin2x 



-C-dyX , 

S+dyX + S-jA^ sin 2x 
7(l + C+sin2x) 



A. 



For a concise notation we introduce the Ginzburg-Landau 
parameter = (A/^o)^ = ^/{^ttKtj^j'^^^) with A = 
[l/(87r7^iir?7Q)]^/^ as the London penetration length (we 
take h = 1). The third equation leads to Ay — > —Ay 
under {9, a) — > {—9,— a). As a consequence we ob- 
serve that the free energy is invariant under the oper- 
ation {9, a, Ay) — > {—9, —a, —Ay), which will be impor- 
tant again later in the analysis of the Josephson effect. 

These equations will now be solved numerically, al- 
though we will also present in |^ an approximate analyt- 
ical solution. For this purpose it is useful to turn to the 
dimensionless variables, measuring lengths in units of ^Oi 
y — y/£,Q, and using for the vector potential = 7Co^j- 
The free energy per unit area of the domain wall can then 
be written as 



Co 



dy 



1 Q ■ 7r. 

- + — sin' 2x 

2 4 ^ 



(dyX)' 



+C+ sin 2x i 



(%X)^) + 2S^axay sin2x 



-|-2%x (S+tty - C-ttx) + K^idyaxY" 



The numerical results for the spatial dependence for the 
order parameter and the vector potential derived from x, 
ay and ax minimizing / are shown in |Bl 

In Fig. [5] we show the domain wall energy per unit 
area as a function of a and 9. We do not use here a as 
a variational parameter but search for the variational lo- 
cal minimum for given a. For this calculations we choose 
the anisotropy parameter to be negative, v = —0.5. Note 
that the result for i' — 0.5 follows immediately from the 
symmetry relation, f(a,9,v) = /(a, ± 7r/4, — :^). For 
9 — 7r/4 we observe two degenerate minima of /(a, 9, v) 
at amiri ~ ±0.467r corresponding to the stable domain 
wall configuration for given 9 and i'. For any other an- 
gle 9 this degeneracy is lifted leading to a stable and 
metastable minimum, which are located at a = and 
a — TT, respectively, if 6* = 0. 

The energy of the stable and metastable domain wall 
configuration depends on the angle 9, as can be seen 
in Fig. H where we plot fmin{9) = f{amin{9),9) for 
V = —0.5 (the solid line marks the stable and the dashed 
line the metastable states). We find that a minimum of 
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FIG. 2: Domain wall energy as a function of the phase dif- 
ference a = (f)+ — <j)- for rotations of the domain wall by the 
angles 6 = Q (solid line), 7r/8 (dashed line) and 7r/4 (dotted 
line). The anisotropy parameter is chosen u = —0.5. 
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this energy occurs at an angle 9 ~ ±0.2l7r away from 
9 — Q. Analogously by symmetry such a minimum is 
found at w ±(7r/4-0.2l7r) = ±0.047r for v = 0.5. From 
this we conclude that there are special orientations for 
the domain wall which are energetically favorable and de- 
pend on the anisotropy properties of the superconductor. 
These special orientations need not to lie along symmetry 
axes or symmetry planes. Note that also the angle with 
respect to the 2:-axis is important in this respect, since 
domain walls parallel to x-y-plane are probably most sta- 
ble. However, here we consider only the case of domain 
walls parallel to the z-axis, as they are generally most 
important for the modification of the Josephson effect. 

As a reference we consider also the case of an isotropic 
Fermi surface {v = 0) which naturally does not show any 
dependence on the angle 9. We find that the most stable 
domain wall state corresponds to the phase a = which 
agrees with the result obtained from a corresponding mi- 
croscopic model calculation based on a quasi-classical ap- 
proach il5|. 
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FIG. 3: Energy of stable (solid line) and metastablc states 
(dashed line) of the domain wall obtained through the mini- 
mization with respect to a: fmin(O) = f{amin,6), as a func- 
tion of the angle of the domain wall and for v = —0.5. 
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FIG. 4: Domain wall energy for an isotropic Fermi surface, 
i.e. v — 0. 



FIG. 5: The stable and metastable values amin as a func- 
tion of the anisotropy of the Fermi surface u, for a domain 
wall parallel to the crystal a;-axis, i.e. 9 = 0. The solid 
lines correspond to stable states while the dashed line denotes 
metastable values ametastabie. = ±7r. 

In Fig. [5]we show the behavior of a„im as a function of 
v for the angle 6 — Q. There are two obvious regions; for 
V < Vc = 0.057 the minimum corresponds to amin = 
and for v > Vc'we find two degenerate values. In addition 
metastable states (indicated as dashed lines) appear at 

Ctmetastable = ±7r for V < -0.12. 



C. Domain wall at the surface 

In view of our later discussion of the Josephson junc- 
tions intersected by domain walls we now consider the 
situation of a domain wall ending at the surface of the 
superconductor and being pinned at a defect somewhere 
in the bulk (for illustration we simplify our model intro- 
ducing a columnar pinning center along the z-axis). We 
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now search for the possible domain waU configurations. 
The domain wall has to compromise between being as 
short as possible and as close as possible to the orien- 
tation (angle 9) minimizing its wall energy whereby, in 
principle, both stable and metastable (local minimum) 
situation can play a role. 

We restrict ourselves here to the case of < which 
is relevant for Sr2Ru04 0, Let us start with the 
situation that the surface normal vector n is assumed 
to be along a crystal main axis in basal plane, say the 
X-axis. 
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case of n = (100). Nevertheless, they represent domain 
wall states which are, in principle, accessible depending 
on the history of the system as there are two local en- 
ergy minima in Fig. [SJ a stable (global minimum) and a 
metastable (local minimum). 
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FIG. 7: Stable branch (solid line) and metastable branch 
(dashed line) of the minimum energy for a domain wall at 
the surface : fmin{0) / cosO , for u = —0.5. 



FIG. 6; Two possible configurations of the domain wall: the 
domain wall may jump between the two positions if it is driven 
by an external magnetic field, as shown below. 

The columnar defect along the z-axis is located at a 
distance xi from the surface (Fig. [B]). The domain wall 
energy per unit length in z-direction for a phase a and 
an angle relative to the x-axis is given by 

Ed{a,e) = xj-^^ . (14) 
cos a 

with f{a,d) defined in Eg. ([T4l) . This is a symmetric 
function under {a, 9) — {—a,— 9) and plotting Ed{a,9) 
we find two stable situations with {amin,9min) and 
i^arnin, —9min)- The two Configurations are depicted in 
Figini We show in Fig. [7] the corresponding domain wall 
energies of the stable and metastable branch given by the 
solid line and the dashed line respectively. The lowest en- 
ergy contributions correspond to the two rather shallow 
minima of the stable branch. 

If the surface normal vector does not lie along a main 
axis (but within the x-y-plane) then the situation be- 
comes more complex. While in the previous case two 
degenerate stable domain wall states were found, under 
general conditions only two local minima of the domain 
wall energy exist which are not degenerate. In Fig. [8] 
we show the domain wall energy for an angle 9 = 0.27r 
of the normal vector n relative to the x-axis. Because 
the domain walls for given orientation possesses phases 
a for stable and metastable energy minima, there are 
two branches of Fig. [S] which correspond to possible do- 
main wall configuration at the surface. These configura- 
tions are neither degenerate nor symmetric unlike for the 
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FIG. 8: Stable branch (solid line) and metastable branch 
(dashed line) of the minimum energy for a domain wall at 
the surface for v — —0.5, when the surface normal vector 
makes an angle 9 — 0.2n with the crystal a;-axis (6 — 9 — 9). 



IV. JOSEPHSON COUPLING BETWEEN s- 
AND p-WAVE SUPERCONDUCTORS 

The symmetry aspects of the Josephson coupling be- 
tween a conventional s-wave and an odd-parity supercon- 
ductor have been discussed many years ago [I3!,ll8|. Con- 
sidering the standard lowest-order Josephson coupling we 
are confronted with the apparent two-fold obstacle that 
there is a mismatch in the spin configuration (singlet ver- 
sus triplet) as well as in parity of the orbital part of 
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the two pair wave functions. It has been shown, how- 
ever, that spin-orbit couphng and the reduced inversion 
symmetry at the interface are sufficient to yield a finite 
coupling. The lowest-order Josephson coupling for an 
interface with normal vector n and the s-wave order pa- 
rameter tp and the p-wave order parameter r) has then 
the form 

Tj = -t{n) I dS {?A*(J7 X n) • i + ip{'q* x n) ■ z} , 

J i 

if the d-vector is parallel to 2-axis as for the chiral p- 
wave state. Here t{n) denotes the coupling strength and 
the integral runs over the interface. Unlike the coupling 
between two singlet-superconductors, which can be es- 
timated through the experimentally determined normal 
state tunneling conductance, in our case the difference in 
the pseudo-spinors yields spin-active (spin-flip) tunnel- 
ing processes which depend on the spin-orbit coupling of 
the two materials. For both Pb and Sr2Ru04 spin-orbit 
coupling is not small, such that it can be expected that 
the Josephson coupling is comparable to an ordinary, al- 
though a reliable estimate is not easy due to the complex 
band structure. 

Using J- J as a boundary term in the Ginzburg-Landau 
equations we can derive the following Josephson current- 
phase relation of the junction: 



J 



We assume n 
to 



2eK, 

— — {ip n ■ Vip — c.c. I 
2 t 

— - [V'*(J7 X n) • i — c.c] 
in 



(16) 



(1,0,0) such that the coupling reduces 



J 



2et 



Aet 



rjy I sin if = Jo sin ip , 
(17) 

with Lp — (j)y — <f>g as the phase difference between the 
order parameters of the p-wave and the s-wave supercon- 
ductor. In this geometry only the Tyy-component of the 
p-wave side order parameter contributes to the Joseph- 
son coupling. This is the order parameter component of 
the pairing state whose nodes point towards the inter- 
face. Only this component can combine with the spin to 
conserve the total angular momentum of the Cooper pair 
in the tunneling through the interface. 



V. INTERFERENCE PATTERN 

We want to calculate the Josephson critical current 
assuming that several domain walls intersect the Joseph- 
son junction from the Sr2Ru04 side. In the following we 
model the effect of the domain wall as a step-like phase 
shift at the point of intersection. This approach is justi- 
fied, if we consider the Josephson coupling as weak such 
that the currents are small and the Josephson penetra- 
tion depth Aj defined below is longer than the extension 



of the junction and even much longer than the width of 
the domain walls 1191. 



A. Junction in a uniform magnetic field 

Analogous to ordinary Josephson junctions we find also 
in junctions between an s- and a chiral p-wave super- 
conductor a relation between the derivative of the phase 
ip = (j}y — (j)s with respect to the coordinate along the 
contact and the magnetic field threading perpendicularly. 
Thus with the assumption that the junction normal vec- 
tor is n = (1,0,0), we examine the variation of the phase 
ip and magnetic field Bz along the y-direction and obtain, 

2tt 

dyf{y) = -^deffBziy) + dy(l){v) , (18) 



where the effective width of the Josephson contact (par- 
allel to n) deff = d + Xs + Xp includes the real width d 
and the London penetration depths on both sides, Xs,p 
(for Pb As ~ 35nTO and for Sr2Ru04 Xp ~ 160nm). In 
contrast to the standard case we add here a contribution 
dy(j) which corresponds to the intrinsic phase variation on 
the p-side which, for example, can be induced by a do- 
main wall intersecting the Josephson contact, as we will 
discuss below. However, also faceting of the interface can 
introduce such contributions. 

The spatial variation of the phase ip{y) obeys the ex- 
tended Sine-Gordon equation. 



1 



(19) 



with Josephson penetration depth Aj = 
{c4>o/87r^ Jode//}^^^. Assuming now that Aj is larger 
than the extension L of the Josephson junction along 
the y-direction, we find that, for a uniform external field 
ip is given approximately by 



27r 

dy'P = -^deffH + dy(, 



ip(y) = ky + (t>{y) + P 



(20) 

where k — 2TrHdef j /^q and f3 is an integration constant. 
The total current is then 



+L/2 

I ^ LzJq I dy sm{ky + (j){y) + l3) 

l-L/2 



(21) 



with Lz the extension of the junction along the z- 
direction. This allows us now to discuss the effect of 
domain wall on the interference pattern in the maximal 
Josephson current obtained by maximizing / with respect 
to (3. 



B. Phase (piy) 

Before addressing the interference pattern it is neces- 
sary to determine the change of the phase 4'{y) as we pass 
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through a domain wall along the interface of the Joseph- 
son junction. This phase enters Eq. ([TS]) through the 
relation 

My - = dyCb , (22) 

on the side of the p-wave superconductor. Consider now 
a single domain wall which reaches the junction interface 
at a 90°-angle, i.e. 9 — 0. The phase shift A(p through 
the domain wall can be determined by integrating along 
the 2/-axis, 

Acb{e = 0) = / dy dycby - — / dyAy 

27r /■+°° 

= ±a{e = 0) - — / dy Ay, (23) 

where Ay has been calculated in Eq. ([T4l) . The plus (mi- 
nus) sign of the phase a reflects the two possible situa- 
tions for the domain wall: moving along the y-axis we 
cross a domain wall from r]- — > 77+ (-1— sign) or vice versa 
( — sign). Since the two situations are related through 
the transformation y — > — j/ it is clear that the phase 
shift A(/) is identical in absolute magnitude for the two 
cases. 

In general, this phase difference depends also on the 
angle of intersection 0, the angle between interface nor- 
mal vector n (in x-y-plane) and domain wall. A rotation 
by an angle 9 involves a phase shift e^*^ for rj±. First we 
consider the case ry_ — 7y+ for which A(j)y = <^-f — 0_ = a, 
if n is parallel to the domain wall. Rotating n by the an- 
gle 9, we obtain A(j)'y = - </>'_ = (^+ + 9)- (0_ -9) = 

a + 29 — a' . Thus the phase shift through the domain 
wall is then given by 

Am = «(^) - -J- / dyAy 
= a'{9)^29 

-— / (cos2M;-|-cos6lsinM;)d2/',(24) 
^0 J —00 

where Ay and A'^ are solutions of Eq.(fT4|) for 6* 7^ 0. 

For an interface with n = (100) the most stable states 
are characterized by the two orientations of the domain 
wall ±9min, according to the previous discussion. By 
symmetry, we know that the free energy is invariant 
under the transformation [a, 9, Ay) — (— a, — 0, — Aj,). 
Changing from 9 mm to —9mm and keeping the same 
boundary conditions (+amm(-6'™„) = -amjn(^'mm)), 
the phase shift simply changes sign 

A0(-0„™O = -A(j){9m,n) ■ (25) 

Similarly, going from the configuration ?7(±oo) — rj^. to 
r7(±oo) = 77^ (keeping the same angle 9min), i-e. under 
a time reversal transformation, only the sign of the phase 
shift is changed 

A(j)+i9m^n) A0_(6',„„) = -A(t)+{9m^n) , (26) 



where we have introduced the notation A(j)+ = A(j){ri_ 
?7_i_) and A(f>_ = A(f){r]_^_ — >■ ri_). We note that this is in 
agreement with the fact that the later transformation 
is equivalent to an inversion of the boundary conditions 
done by the transformation a — > —a, which causes a sign 
change for the phase difference. 

Therefore, to each domain wall intersecting the 
Josephson junction we can associate a phase difference 
A(l3+{9min) (resp. A(j)^{9mm)) corresponding to the ge- 
ometry ri{y — ±00) = rjj. (resp. ri{y — ±00) = r]^) 
which depends on the orientation of the domain wall. 

For orientation of the interface normal vector different 
from n = (100), we compare the two lowest energy con- 
figurations of the domain wall (stable and metastable) 
state. Using the above calculation scheme we can show 
that the phase differences are in general different for the 
two states, yielding Acpi and A(j)2- Under time reversal 
operation, changing between the two cases 77+ rj and 
r7_ — >■ r/_|_ , we obtain the opposite sign of the phase shifts 
again. Note that, in general, the two phase shifts, A(/)i 
and A(j32 are comparable in magnitude with those cases 
where n lies along a symmetry axis. Note, however, that 
there is an energy difference between the two domain wall 
states, if they are not degenerate. This energy difference 
depends on the depth over which the domain wall is in- 
fluenced by the change of configuration at the interface. 



C. Model for intersecting domain walls 

We extend now our model to a Josephson junction 
which is intersected by several domain walls. First con- 
sider the case of one domain wall reaching the interface 
at the position y = yi, which yields the Josephson phase, 

^(y) = ky + (j){y) + p = ky + A(j)i^Q{y - yi) + p . (27) 

Here A(f>if^ denotes the phase shift at the domain wall 
as calculated in the previous section. The index /i labels 
the two types of domain wall state (stable or metastable 
minima). It is admissible to use a step function 0{y) to 
describe the spatial change of the phase, since the exten- 
sion of the domain wall is small compared to the length 
of the Josephson junction and the Josephson penetration 
depth Aj. This leads to a piecewise constant phase shift. 

The generalization to an array of domain walls is 
straightforward. For N successive domain walls, the 
Josephson phase difference across the junction is given 
by 

N 

^{y) = P + Y,A(j>,^^ e{y-y,) + ky , (28) 

where </'(?/) depends on the configuration 
{/ii, /i2, ■ ■ • , Miv} of all domain walls. We will now 
show that the distortion of the Josephson interference 
pattern through ip{y) depends on these configurations. 
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D. Modified interference pattern 

The total Josephson current which traverses the junc- 
tion is given by Ea. (pT|) . Using our piece- wise constant 
approximation for (/^{y) we obtain readily the following 
expression, 



N 



ffi = ^ Vsin /? + Va</>,,. + !L£(^^±i±M 

'7r$ (y,+i - y. 



$0 



L 



(29) 



with A(t>Qf,g = 0, yN+i = L/2, yo = -L/2 and Iq = SJq 
{S = LzL being the interface area), and kL^Q = 27r$. 

The Josephson critical current, i.e. the maximal su- 
percurrent possible for this configuration, is obtained by 
the maximization of / with respect to /3. This condition 
reads 



5/ 

dp 



leading to 



with 



N 



arctan 



(30) 



(31) 



N 



sm 




X sm 



_ (y»+i + Vi ) 
$0 L 

(y^+i - y^) 



$0 L 

The Josephson current can be written as 

— = -^{Ai sin/3 + A2 cos/S} , 



(32) 



(33) 



and finally we find 



7r<i> 

$0 



{Ai Sin^,„ax + A2 cos /3,nax) 



7r<i> 



Al + Al 



(34) 



In Figini and FiglTO] we show two simulations for in- 
terference pattern of the Josephson critical current in a 
magnetic field, using two random domain wall configu- 
rations (indicated in the inserted panels) where we as- 
sumed A'^ = 10 for the number of intersecting domain 



walls. The randomness occurs in the positions of domain 
walls as well as in the sequence of {/ii,...,/xjv}. The 
deviation from the standard Fraunhofer pattern is obvi- 
ous. However, the maximum of the Josephson current 
lies rather close to H — 0. This feature has to do with 
the fact that, in the case of generic random configura- 
tions, the phase shift does not vary much overall, i.e. 
10(2/ = L/2) — 0(y = —L/2)\ <^ 2tt. Stronger deviations 
can be observed, if we bias the domain wall configuration 
in a way as to have a larger net shift, e.g. by assuming 
for all domain walls the same sign of A0. 



= 0.6 - 



0.4 




FIG. 9: Josephson critical current as a function of the exter- 
nal magnetic flux for ten intersecting domain walls randomly 
configured with a total phase shift 2|A</!>|. The anisotropy 
parameter is chosen u — —0.6 fl6l. [l7j. 



= 0.6 - 




FIG. 10: Josephson critical current as a function of the exter- 
nal magnetic flux for ten intersecting domain walls randomly 
configured with a zero total phase shift. The anisotropy pa- 
rameter is chosen ly = —0.6 [l^. [l7j. 



E. Hysteresis and noise effects 

An important feature supporting the idea that domain 
walls are involved in producing the irregular interference 
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pattern is the observation of hysteresis effects by Kid- 
wingira et al, when the external magnetic field was cycled 
between positive and negative maximal fields [l|. They 
argued that the applied magnetic field induces a rear- 
rangement of the domain walls in the sample, and they 
substantiated their claim by simulations looking at the 
effect of shifted domain wall positions. Since the domain 
walls are pinned at defects of the sample large shifts in 
positions are rather unlikely. 

Hysteresis effects are rather easily discussed within our 
simplified model in Eq.(I2Hl HSl IM])- The free energy of 
the junction in a magnetic field can be approximated by 

F (/3,$,{^i,...,/iA,}) 

(35) 

We neglect the change of the intersection points yi for 
different configurations and examine the condition to 
minimize the free energy. Moreover we assume that the 
interface is smooth and that we have the same two es- 
sentially degenerate configurations for all domain walls: 
A</)a = A(/)+ > and A(/)i2 = A0_ < 0. Thus, n 
is assumed to lie very close to a high-symmetry axis of 
Sr2Ru04, say n = (100). For values $ > the free 
energy in Eg. ([55)1 can be lowered by choosing A0i^. = 
A(j)-, because under this condition 



striking feature is the shift of the maximum of the crit- 
ical current. These curves are obtained assuming that 
changes of domain wall configurations during the field 
sweep, when the critical currents is measured, can be 
neglected. Reorganizations of the domain wall config- 
urations during the field sweep would most likely lead 
to discontinuities in the critical current. This kind of 
behavior is also observed in a set of measurements by 
Kidwingira et al (Fig.3B in [ij) where the field sweeping 
range is restricted to rather small fields only. Also in this 
case a hysteretic behavior of critical current (dependence 
on the field sweep direction) was observed, though less 
pronounced. It is not unlikely that some part of the noise 
on these data can be interpreted as an effect due to the 
domain walls. 
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(36) 



Note that there are special values of $ where the full "po- 
larization" of A^ may not be the best choice. However, 
considering a sweep of the field to maximal value ^max 
and back would favor Ea. p6|) . The same argument can 
be used for negative fluxes, driving the domain walls to 
adopt A(j)ifj_. = A0_|_. Note that for metastable domain 
wall states a "polarization" would only be possible, if the 
energy expense of the metastable configurations is suffi- 
ciently small. We consider now geometries for which this 
is true. 

Starting at zero field, the domain walls shall be essen- 
tially randomly configured. Such a situation would lead 
to a critical current pattern as shown in Fig. [TOl If the 
magnetic field is increased up to a sufficiently high value, 
then the domain walls would likely polarize after a cer- 
tain waiting time. When we decrease now the field to 
measure the interference pattern, we observe the inter- 
ference pattern modified by the polarized domain walls. 
After reaching a sufficiently large negative field value, the 
domain walls polarize in the opposite way. Therefore, 
tuning the field back towards positive values we find an 
altered interference pattern of the critical current. Simu- 
lation results taking the two (polarized) domain wall con- 
figurations into account are shown in Fig llll The most 



FIG. 11: Shift of the maximum of the critical Josephson cur- 
rent for positive field (solid line) and negative field (dashed 
line) as observed by Kidwingira aZ [l[|. 

Turning back to the hysteresis effect with polarized do- 
main wall configurations, we might use the shift of the 
position of the maximal critical current in order to es- 
timate the density of domain walls. We denote the fiux 
of the maximal current as which can be determined 
approximately by the condition 



N 

E 

i=l 



A0,,,, = -(0(L/2)-(/.(-L/2)) = -NAcI) 

(37) 

assuming for the last equality a completely polarized 
domain wall configuration where all domain walls con- 
tribute the same phase shift Acf). First, we conclude that 
the maximum lies at $mc > for a sweep down from 
the positive field side (opposite for the opposite sweep 
direction) in accordance with experimental findings [l|. 
The number of domain walls intersecting the Josephson 
junction is then given by 



N : 



27r$, 



A0$o 



(38) 



Kidwingira et al find in a measurement of hysteresis effect 
the magnetic field Bmc « 0.8G (Fig.3A in [3j) which 
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deffLBmax » 16 X G^m 



0.8$o- 

100/im, As = 



yields a flux 

For this estimate of we took L 
35nm, d = lOnm and Ap = 160nm. 

Using Eq. with = —0.6 wc find A(/) in the range 
of 0.2 - 0.5 depending on the orientation of the domain 
wall with respect to the interface. From this we obtain 
= 10 - 25 which leads to a mean distance between 
domain walls of ~ 4 - lO/im. 

Kidwingira et al also report the presence of peculiar 
noise in the time dependence of the voltage drop of the 
junction in a constant current slightly above the critical 
current [ij. The question arises whether this feature can 
also be attributed to the dynamics of domain walls. As 
we have seen above, the modification of the domain wall 
configuration {/ii, . . . , /iA?} changes the critical current. 
The current-voltage characteristics of a Josephson junc- 
tion for currents / immediately above the critical current 
Ic is very non-linear. The standard textbook form of 
the current-voltage relation of a Josephson junction gives 
usually a good approximation of the general behavior 



V = Ry/P - /2 



(39) 



for I > Ic- For a given current / the time dependence of 
the critical current, Ic{t) = Ic + 5Ic{t) leads to 



6V{t) w Rh 



Slcjt) 



(40) 



which can be large for I close to Ic- The reported noise 
feature (Fig.SD of Q) suggest that the critical current 
fluctuates between two values, which could correspond 
to two domain wall configurations (see. Fig fT2)) . 



the s-wave superconductor and Sr2Ru04. The effect of 
faceting is two-fold. First, a faceted surface can lead 
to additional stable and metastable domain wall config- 
urations as domain walls can be pinned at surface in- 
homogeneities. Thus, they add to the number of pos- 
sible junction states. Second, faceting corresponds to a 
varying interface normal vector. It is easy to see that 
the modulation of the inplane normal vector angle 6 by 
59{jj) corresponds directly to the Josephson phase varia- 
tion 5(t){y) = 89{y) {\59\ <C tt), if the length scale of the 
faceting is larger than the coherence length of Sr2Ru04 
(^ ^ 80nm). As a random phase modulation of this 
kind with \5(j){L/2) - 5(j){-L/2)\ < Wl does not lead to 
a significant shift of the maximal critical current in the 
interference pattern and no hysteretic effect is possible, 
faceting alone cannot be responsible for the features ob- 
served. Together with a spatial variation of the Joseph- 
son coupling the variation of the phase due to faceting 
can contribute to the static random structure of the in- 
terference pattern. Fine structures in the interference 
pattern (variations of the critical current on small field 
changes) are usually due to variation on long length scales 
along the junction and are also known for conventional 
Josephson contacts. Note that the effect of faceting in 
d-wave high-temperature superconductors has a stronger 
impact on the interference pattern, since the phase jumps 
between and tt can lead to phase wandering along the 
junction such that (/)(L/2) — (j){—L/2) > 2tt and produces 
a field shift for the maximal Josephson current, as ob- 
served in experiment [ll|. 



V 




VI. LONG JOSEPHSON JUNCTION SOLUTION 

For completeness we address here a further interesting 
feature connected with a domain wall intersecting the 
interface. This time, however, we change to limit of Aj 
being much shorter than the extension of the junction. 
In order to describe the behavior of the phase we have 
then to solve the Sine-Gordon equation 



a2(p = -^sinv5 + 9>. (41) 



FIG. 12: Schematicl-V characteristic of a Josephson junction. 
Two values of the critical current are assumed. Id and /c2, 
which for the same applied current lead to different values 
of the voltage on the junction, Vi and V2, respectively. A 
time sequence of changing between the two critical current 
values (by change of domain wall configurations) gives rise to 
steplike noise on the voltage signal. 

Our theoretical discussion of the interference pattern 
did not include the faceting of the interface between 



The solution here has a kink shape which is slightly mod- 
ified by the last term of the right-hand side. This term 
acts like a source term in the vicinity of the domain wall. 
Without explicitly writing down the solution it is clear 
that within the length scale of Aj from the domain wall 
ip{y) approaches a constant value. As discussed earlier 
the derivative dyip corresponds to a local magnetic field, 
such that here a well-localized magnetic flux line appears 
at the position where the domain wall meets the interface. 
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The enclosed magnetic flux can be easily determined by 

rL/2 



[Ay{+a,y) - Ay{-a,y)]dy 

-L/2 



L/2 



27r 



{dyip{y) ~ dy(j){y)) dy 



$0 

2tt' 



n27r 



$0 

27r 



(42) 



where n is a positive integer. We finally see that the 
minimal possible fluxes is given by 



"2^ 



A^, 
"27- 



$0 



(43) 



Since ^ < 1, we encounter here fractional vortices which 



are generally a sign of broken time reversal symmetry 
[3, H^l ■ The detection of such well- localized fluxes could 
be used to detect the position of domain walls. 



VII. CONCLUSION 

Motivated by the experiments of Kidwingira et al. on 
Josephson junctions between Sr2Ru04 and the conven- 
tional superconductor Pb [l| , we studied effect of domain 
walls on the Josephson interference effect in a magnetic 
field assuming that Sr2Ru04 is a chiral p-wave supercon- 
ductor. For this purpose we analyzed the domain wall 
structure and showed that its internal phase structure is 
crucial for the Josephson effect, if domain walls intersect 
the Josephson junction. The anisotropy of the electronic 
band structure plays an important role, as it influences 
the phase shift of the superconducting order parameter 
between the two types of chiral domains. In addition, it 
determines the energetically most favorable orientation 
of the domain wall. 

The theory presented here is able to explain the most 
important features reported by Kidwingira et al [l|. 
(1) The interference pattern in a single junction devi- 
ates from the the standard Fraunhofer pattern not only 
through irregularities, but also shows a distinct asym- 
metry between positive and negative magnetic fields and 
the maximum of the critical current can be shifted away 
from zero field. This can be attributed to the phase 
shifts introduced by domain walls intersecting the inter- 
face between Sr2Ru04 and the conventional supercon- 
ductor. (2) Kidwingira et al report strong differences 
between the interference pattern of different samples. 
This may be explained by the fact that the samples were 
prepared in different ways, with possibly also different 
normal vector directions as well as different degrees of 
faceting. (3) Cycling the field continuously covering a 
positive and negative field range, a hysteretic behavior in 
the interference pattern appears which can be understood 
as a field-driven motion of the domain walls, whereby 
the anisotropy of the domain wall energy could play an 
important role, as our model shows. Allowing for the 



change of domain wall configurations in time we can also 
understand the noise effect seen in the current voltage 
measurements for currents above the critical current. 

The question arises whether it would be possible to 
test the domain wall scenario by generating a single- 
domain superconducting phase in Sr2Ru04. This would 
be most straightforwardly realized by field-cooling. Un- 
fortunately, this procedure might lead to vortex trapping 
in both superconductors which would jeopardize a clear 
outcome. Domains naturally arise in the zero-field cool- 
ing process, as superconductivity nucleates in the sam- 
ple with some extent of inhomogeneity, since the chiral 
p-wave state is rather susceptible to disorder effects. In 
this way both chiral states can emerge in the sample in 
different regions and the domain walls appearing between 
them are eventually pinned in the sample, making it dif- 
ficult to anneal the sample to a domain wall free phase. 
It is difficult to estimate a priori the density of domain 
walls. However, our analysis suggests that domain walls 
intersecting the Josephson junctions in the experiment 
may be separated by a few micrometers which is con- 
siderably larger than the extension of the domain walls. 
This estimate may be interesting in the context of recent 
measurements by Kirtley et al aiming at the observation 
of spontaneous magnetic flux at surfaces of Sr2Ru04, as 
expected for chiral p-wave states P]. Since these exper- 
iments gave a negative result for the presence of such 
fluxes, it was speculated that the dense population of 
chiral domains would cancel out the signal for the scan- 
ning SQUID microscope used. Our estimate of domain 
wall density corresponds roughly to the spatial resolu- 
tion of the SQUID microscope. Thus it remains unclear 
whether domain walls really could explain the absence of 
a positive signal from spontaneous flux. Extensive theo- 
retical studies on this issue point to an important puzzle 
in this context ^l|. 

Eventually, we would like to note that this experi- 
ment does not decide between even- or odd-parity pair- 
ing. While the experimental situation satisfies the se- 
lection rules for the coupling between the chiral p-wave 
and an s-wave superconductor, the effects observed in the 
Josephson interference experiments would be similar for 
a chiral d-wave superconductor such as the kz{kx ±iky)- 
wave state belonging to iJ^-representation of the tetrago- 
nal point group -D4/1, for which, however, additional com- 
plications appear due to the nodal gap structure. At the 
present stage, however, we do not have any experimental 
signature in Ref. [l| , which would rule out the E'g-state. 
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Appendix A: Analytical approximation 



We start from the free energy functional (|T3)) where we 
neglect the terms in A^- Varying with respect to x 
Ay gives, respectively, 



-.m4x~ S+-f{dyAy) 



+ C+[cos2xh'Al + (dyxf) + sm2x{d'yx)m) 

_ S+dyX 

- 7(l + C+sin2x)- ^^'^ 

Following the argument of Ref . [l3| we neglect the con- 
tribution of the third term of (|A1[) and the sin 2% term 
in the denominator of Ay in (IA2p . This leads to the sim- 
plified equation for x{y)y 



(A3) 



where Q — Q/[£,q{1 — S^)]. This has the following kink 
solution 



= arctan c 



-VQi 



(A4) 



which we use to determine Ay. Inserting x and Ay into 
the free energy functional we derive an analytical form 
of the variational domain wall energy per unit area. 



f{a,9) = Kr^lJQll- 



sl 



Si 



C.TT 



arctan 



CI 



1-Cl 



1 + C4 



k5) 



In Fig. [13] we show the domain wall energy as a func- 
tion of a for the different angles 9 — 0, tt/8, 7r/4, with 
the anisotropy parameter chosen u = —0.5. We see that 
those results are qualitatively and even quantitatively 
close to the numerical results plotted in Figl2] In Fig. 
[14] we plot the free energy density at = 0. Here some 
short-comings of the above approximation becomes obvi- 
ous, since there a two minima of the energy in contrast 
to the single one of the numerical result shown in Fig. 
[3] This discrepancy originates from neglecting the vec- 
tor potential A^ and the corresponding terms in the free 
energy. 



Appendix B: Numeral solutions for u 



-0.6 



We show here the numerical solutions of the structure 
of the domain wall in the most stable state for i' = —0.6 



and compare with the analytical approximation. We plot 
in Fig. [Tn]the function xiu) where the solid line is the 
numerical solution and the dot-dashed line the approx- 
imative solution in Eq. (|A4p . The deviation is rather 
1.6 
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FIG. 13: Domain wall energy (derived from the analytical 
approximation) for a chosen anisotropy parameter 1/ = —0.5 
as a function of the phase difference a = — 4>- when the 
angle of the domain wall with respect to the crystal x-axis 
are: 6* = 0, 7r/8 and 7r/4. 




FIG. 14: Domain wall energy derived from the analytical ap- 
proximation for an isotropic Fermi surface. 



small. In Fig. [TB] we show the two components of the 
vector potential, the solid line for the numerical and the 
dot-dashed line for the approximate solution. Here the 
discrepancy is larger. 

Fig. [T7] and [15] depict the modulus of both order pa- 
rameter components, \'qx\ and \i]y \ and the phase of the 
order parameter passing through the domain wall (we 
choose (f>± = ±a/2). Obviously, the two phase differ- 



ences follow the relation A^^; 



TT and A(j)y = a. 
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FIG. 15: Spacial dependence of xiv) through a domain wall, 
calculated numerically (solid line) and with the analytical ap- 
proximation (dot-dashed line) for i/ = —0.6. 
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FIG. 16: Spacial dependence of the vector potential through 
a domain wall. Numerical solutions of Ay (solid line) and A^c 
(dashed line). Analytical solution of Ay (dot-dashed line). 
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FIG. 17: Spacial dependence of the modulus of the order 
parameter through a domain wall : \rix\ (dashed line), \riy\ 
(solid line). 
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FIG. 18: Spacial dependence of the phase of the order pa- 
rameter through a domain wall : 4>x (dashed line), cjiy (solid 
line). 
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